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We study the transport through evanescent waves in graphene quantum dots of different geome- 
tries. The transmission is suppressed when the leads are attached to edges of the same majority 
sublattice. Otherwise, the transmission depends exponentially on the distance between leads in 
rectangular dots, and as a power law in circular dots. The transmission through junctions where 
the transmitted and reflected currents belong to the opposite valley as the incoming one depends 
on details of the lattice structure at distances comparable to the atomic spacing. 
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I. INTRODUCTION 

Since its synthesis and identificationiii^, graphene has 
attracted a great deal of interest, because of its unusual 
fundamental properties, and possible applications^iii^. 
Being purely planar material (actually, graphene is 
the first example of truly two-dimensional crystals) 
and demonstrating a very high charge carrier mobility 
graphene is considered as a perspective base for a post- 
silicon electronics. 

The carriers in graphene are described by the mass- 
less Dirac equation and possess a pseudospin degree of 
freedom (which is, actually, a sublattice label) and chi- 
rality related with it. A number of unusual properties 
follow from this, and, in particular, the existence of local- 
ized states at the Dirac energy (that is, at the neutrality 
point), and a new transport regime dominated by evanes- 
cent waves, also at the Dirac energj*^ (see also the ex- 
perimental studies in^'^). Midgap states with the energy 
close to zero (further we will count the energy from the 
Dirac energy) were initially found at graphene edges with 
a perfect termination in one of the two sublattices which 
define the honeycomb lattice^ (zigzag edges), and later 
they were generalized to other defects, such as cracks^^, 
vacancies^, generic surfaces with a majority of atoms of 
one sublatticei^, and ripples^. These midgap states have 
similar wave functions to the evanescent waves which me- 
diate the transport in clean graphene when there are no 
charge carriers and the chemical potential coincides with 
the neutrality point^i. The combination of evanescent 
waves and localized states can even enhance the conduc- 
tivity in graphene with defects^ii^. It can also be ex- 
pected that midgap states dominate the transport prop- 
erties of graphene quantum dots^^ which are the subject 
of intensive study now^ '^^i^^'^°i^^ . 

In the following, we extend the analysis the trans- 
port properties of clean graphene with chemical potential 
equal to zero^ to graphene quantum dots of different 
geometries. We first present the model, and then we 



consider three different cases, namely, rectangular dot, 
circular dot, and the corner between two facets represen- 
tative of broad classes of quantum dots. In particular, 
we will demonstrate that for the case of circular dot the 
conductance is very sensitive to magnetic field which may 
give an insight for development a new type of magnetic 
sensors. 

The main conclusions can be found in the last section. 



II. THE MODEL 

We will consider ballistic quantum dots etched from a 
single layer graphene flake, connected by graphene leads. 
We assume that the chemical potential in the external 
reservoirs and the leads is far from the Dirac energy, while 
a gate fixes the chemical potential of the quantum dot at 
the Dirac energy, e = 0. As there are no extended states 
at e = in the quantum dot, transport processes take 
place through evanescent waves induced by the contacts. 
This combination of chemical potentials is similar to that 
considered in^ii for the case of bulk graphene. 

We assume that the leads have a width I. We also 
assume in the following that the carrier concentration 
in the dots is higbj!§.ii, so that kpj^l :s> 1, where kpi^ 
is the Fermi wavelength in the leads. The transmission 
through the constriction induces changes in the momen- 
tum parallel to the interface, ky, of order Aky '^ l~^. 
As discussed below, the transmission through the dot 
is determined by evanescent waves with decay length 
'^^^ <^ ^- The allowed parallel momentum of these modes 
is ky r^ K r^ l~^. Hence, the relevant incoming modes 
in the leads also must have ky = kF]^sm{9) « l~^. The 
relevant states in the leads are focused in the forward di- 
rection, 9 < l/{kpi^l). Within these approximations, the 
incoming and outgoing waves can be written as: 
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where x is the coordinate along the axis of the lead, the 
index i = in, out defines the incoming and outgoing leads, 
A and B label the two sublattices in graphene, and the 
two signs correspond to waves moving in opposite direc- 
tions. The states in cq.([T|) are the same as those used 
in^ii. In the following, we will neglect inter-valley scat- 
tering in the dot, so that the transmission through the 
dot can be analyzed for each valley separately. 

The wave function in the incoming and outgoing leads, 
near the contacts to the dot, Xj = 0, can be written as: 
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where R and T are the refiection and transmission am- 
plitudes, and the carriers in the leads are supposed to be 
electrons (e > 0). 

We assume that surface states exist at the edges of the 
dot. This is the general case, as these states always exist 
when the two sublattices of the honeycomb structure are 
not compensated at the edgei^. The existence of these 
states leads to a depletion of the density of extended 
states for energies |e| <C v-p/L, where L is a typical di- 
mension of the dot. Hence, within this range of energies 
we need only consider the midgap states within the dot. 
Defining the complex variable z = x + iy where x and y 
are Cartesian coordinates, the wave function of the states 
with zero energy must be of the form.®: 



"^K'ix^y) 






K 



{x,y) 






Ik{z) 
9k{z) 

Ik'{z 
9k' {z 



(3) 



where fK,gK, fx' and gK' are analytical functions, and 
K and K' label the two inequivalent valleys in graphene. 
In the following, we will consider only one valley. The 
results can be extended to the other valley in a straight- 
forward way. For the study of the two cases shown in 
Fig. [1] , we neglect processes which lead to the hybridiza- 
tion of the two valleys. Note that we are considering 
ballistic systems, where intervalley scattering can be ne- 
glected. A different situation is that of a 60° edge, where 
intervalley scattering is crucial, even in the ballistic limit. 



FIG. 1: (Color online). Sketch of the geometries of a graphene 
quantum dot analyzed in the paper, a) Rectangular dot. b) 
Ring shaped dot. c) 60° edge. 



For simplicity, we assume periodic boundary condi- 
tions along the vertical direction, as in Ref.®. For a 
graphene ribbon of width W ^ a, the solutions for 
open boundary conditions can be written as superposi- 
tions of waves with opposite momenta which are pos- 
sible solutions of an equivalent problem with periodic 
boundary conditions. The difference between the two 
choices for a ribbon with TV unit cells of length a lies in 
the set of allowed momenta, k± = mT/[{N + l)a],n = 
1, • • • ,N for periodic boundary conditions, and k±_ — 
zL2Tm/{Na), n = 0, ■ • ■ , N/2. If the transmissions to be 
calculated are a smooth function of the momentum k± , 
the results obtained using either choice of boundary con- 
ditions converge to the same value for A^ :^ 1. 



For the rectangular dot considered here, the x coor- 
dinate lies within the range —L/2 < x < L/2, and the 
y coordinate can be written &.& y = W x 6'/(27r), with 
< 9 < 27T. The wave functions inside the dot can be 
expanded using the basis: 
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III. RESULTS 



A. Rectangular dot 

We first consider a rectangular dot, determined by its 
length L, and width W, as sketched in Fig.[l^]. Each of 
the leads be connected to an edge terminated in the A 
sublattice, and the other be connected to an edge termi- 
nated in the B sublattice. 



The incoming lead is attached at the position zi = —L/2, 
and the outgoing lead is attached to the position Z2 = 
L/2 + iW X 9o/{2n). The contact averages out the details 
of the wave functions in the dot over a length of order I. 
Hence, the description of the effects of the contact does 
not require the infinite sum in Eq.(|3]). In the following, 
we include an upper cutoff, rimax ~ W/l. The boundary 



conditions at zi and Z2 ai"^: 
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The zigzag boundary conditions at points at the edge 
with X = —L/2 other than the position of the lead, 6 = 0, 
and at the edge with x = L/2 for 9 ^ Oq read: 

n — — njnax 
n^nax 

Y a„e^"-^/^e™'' =0 61 7^ 6*0 (6) 

An Ansatz which is compatible with the two sets of 
boundary conditions is 
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which leads to the equations 
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and 



r = 



induces a change in the allowed values of the transverse 
momentum, ky, with respect to the case of zigzag bound- 
ary conditions. The momentum shift is^ Aky ^ W~~^. 
An incoming valley polarized mode can be written as a 
superposition of exact eigenmodes in the presence of val- 
ley mixing at the boundaries, with a spread in momentum 
of order Aky ^ l~^. 

As discussed above, the transmission through the dot is 
determined by evanescent waves with decay length k~^ ~ 
I. The parallel momentum associated to these states is 
ky ^ K ^ l^^. This superposition of valley polarized 
states will be changed by the existence of valley mixing 
at the top and bottom boundaries. The effects are of 
order e^'^'^^ ^ e"^/*^, which are exponentially smaller 
than the transmission coefficient obtained above, \T\ ~ 

g-47rL/i^ 

The previous analysis can be extended to a rectangular 
graphene quantum dot in a constant magnetic field. The 
wave functions in Eq.Q in this case become 
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(9) 
The sums in Eq.® are dominated by the terms with 
n ~ irimax- Hence, we obtain: 
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Finally, when the two contacts are attached to the same 
type of zigzag edge, the boundary conditions lead to ei- 
ther a„ = (contacts attached to an A terminated edge) 
or bn = (for a B terminated edge), and, as a result, 
i? = 1 and T = 0. 

The assumption of a rectangular shape implies that 
the top and bottom boundaries may be of the armchair 
type, leading to intervalley mixing-S. These boundaries 
are separated by a length W ^ I. The valley mixing 
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where Ib = \/hc\e\B is the magnetic length, B is the 
magnetic induction. The same manipulations described 
above allow us to calculate the transmission coefficient, 
which turns out to be unchanged with respect to Eq. PU)) . 
The lack of dependence of the transmission on the ap- 
plied field is in agreement with the insensitivity of the 
bulk transport to the magnetic field when the chemical 
potential is at the Dirac point^. 



B. Circular dot 

We now consider a ring-shaped dot, with the contacts 
attached to the inner and outer edges, which are as- 
sumed to be A and B terminated, respectiveljsiZ,. As 
discussed there, a graphene dot with an approximate 
'1 circular shape can be bounded by different edges where 
the majority sublattice is of either type. The advantage 
of assuming the same boundary conditions throughout 
the edge is that the problem retains circular symmetry, 
and the solutions have a well defined angular momen- 
tum, I, where I = i,---lmax ~ R/a, where R is the 
radius of the dot and a is the lattice constant. For 
R/a ^ 1 , superpositions of wavefunctions with I 3> 1 can 
be built, with a small spread in angles, ^(r, 0) j^ only if 
9()-Ae<0<Oo + A9, with A9 < 27r. These states are 
solutions of the Dirac equation at zero energy which are 
not affected by the global properties of the boundaries 
of the dot. Hence, they can be used to describe trans- 
port between positions at the dot edges which are close 
enough so as to be insensitive to the global properties of 
the edges. The case when the transport properties are 



influenced by the type of edges requires input about fea- 
tures at distances comparable to the lattice spacing, and 
will be considered in the next subsection. 

The outer and inner radii are i?i and R2, as shown 
in Fig.[T]3]. The modes at e = inside the dot which 
satisfy the boundary conditions can be obtained by a 
conformal transformation which changes the rectangular 
(cylindrical) dot considered in the previous subsection 
into a ring. This transformation is 
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and e^'^^/^ — R1/R2. Using this transformation, the 
wave functions in Eq.(|3]) become 
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The size of the contact, I, induces a maximum value n, 
nmax ~ 2TrR2/l. From the wave functions in Eq. (ll3p . and 
using the same analysis as for the rectangular dot in the 
previous subsection, we wind, in analogy with Eq. (|10p : 

If we write 2ttR2 — W and i?i « i?2 + L, this expression 
reduces to Eq.HU]) when Ri r^ R2 ~> Ri - R2- 

As in the case of the rectangular dot, the only allowed 
solution when the two contacts are attached to the same 
boundary is i? = 1 and T = 0. 

Dots of different shapes can be obtained using other 
conformal transformations. For instance, the circular dot 
analyzed here can be turned into a dot with a corrugated 
boundary using the mapping 



w (w) — w + Xw^' 



(15) 



where m fixes the period of the corrugation, and A ex 
6R/R™ gives its amplitude. The transmission coefficient 
in eq. p^ acquires corrections of order {R/l){SR/R)"^. 

This analysis can also be extended to a finite magnetic 
field. In that case the wave functions are 
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The boundary conditions imply 
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possible solution is of the form 
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leading to the equations 
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Thus, the transmission coefficient is 
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where $ is the magnetic fiux per ring and $0 = 7r?ic/|e| 
is the flux quantum. The dependence of transport 
through the quantum dot on the magnetic field is ex- 
ponential, rather than oscillatory as in case of the stan- 
dard Aharonov-Bohm effect^^, which reflects a specifics 
of transport via evanescent waves. 

Thus, the conductance of circular quantum dot turns 
out to be very sensitive to the flux through the ring. 
This can be potentially interesting in light of develop- 
ment of magnetic sensors for measurements of very low 
fields without use of superconductivity. 



C. Corner betAveen two facets 

We finally consider a 60° angle boundary between an 
edge with A termination, and an edge with B termina- 
tion, as sketched in Fig.[TJ:]. The leads are attached at 
distances li and I2 from the corner. The transmission 
along a ribbon with this geometry has been analyzed 
numerically in Refsi^i2£iS2. A full solution cannot be 
obtained within the continuum approximation, as in the 
previous cases, since both transmission and reflection at 
the wedge require Umklapp processes changing the val- 
ley index of the incoming electron. In this respect, the 
problem is similar to that of the transmission across a 
p-n junction near the Dirac point^'*. 

We consider the ribbon with an angle shown in Fig. [2] . 
The system has reflection symmetry around an axis 
which connects the vortices of the angles at both sides. 
Far away from the corner, the system reduces to a zigzag 
ribbon. We assume that the edges of the ribbon are at 
y = ±W/2, and W ^ a, where a is the lattice constant. 
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FIG. 2: (Color online). Sketch of the 60° junction analyzed 
in the text. Top left: an incoming wave from the left, built up 
from eigenstates in the vicinity of the K valley has to be either 
transmitted or reflected as a superposition of states from the 
K' valley. Bottom: Scheme used to analyze the transmission 
in the honeycomb lattice. The junction has a reflection sym- 
metry around its center, shown as dashed line. On each side, 
the nearest neighbor hopping problem can be written as a 
sum of coupling between nearest neighbor transverse stripes. 
See text for details. 
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FIG. 3: (Color online). Reflection at a 60° junction like the 
one sketched in Fig. [2] for different junction widths. A^ gives 
the number of unit cells across the junction. The width of 
each ribbon which make up the junction is 3/2aN, where a 
is the length of the C-C bond. The energy is in units of the 
nearest neighbor hopping. 



The wave function of an incoming wave can be written 
as 



*5?(x,y)^ 
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where k and the energy of the state are given by: 
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(23) 

The solution dUD is vaUd for k > W'^. 

The reflected and transmitted waves can be written 
as in Eg. ([2T|) . except that, for the reflected wave the mo- 
mentum is reversed, and the wave packet is built up from 
electrons at the K' valley: 
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(24) 
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FIG. 4: (Color online). Weights of the wavefunction passing 
through the junction sketched in Fig. [2] on evanescent waves 
of different decay lengths. The width of the junction is 16 unit 
cells, and there are 30 evanescent waves. The decay length is 
in units of the lattice spacing in the direction parallel to the 
edges (V3 x a, where a is the C-C bond). Diamonds show 
waves symmetric around the junction, and hexagons show 
antisymmetric waves (note that the waves studied here decay 
j^^as function of the distance to the junction, not to the edges 
of the ribbon). 



This wave function has to be matched to that in Eq. (|2T|) 
at the interface. In the continuum hmit, the overlap be- 
tween wave functions at K and K' valleys is zero, and 
the matching cannot be carried out. We analyze the cor- 
rections induced by the finite lattice by dividing each 
ribbon reaching the junction into stripes, as sketched in 
Fig. [2]. Using a nearest-neighbor tight-binding model, 
each one-dimensional ribbon is coupled to its two near- 
est neighbors, and the Hamiltonian can be written as a 
sum of iV X TV terms, Ti and T2, where N is the number 
of atoms within each ribbon: 
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where M is a stripe index. 

Using the inversion symmetry of the junction, we can 
analyze separately states which are even and odd with 
respect to the axis which joins the two angles of the junc- 
tion. Within each subsector, the problem is reduced to 
the reflection of an incoming wave of energy et at a sharp 
boundary at the location of the junction. The coupling 
between the stripes at each side of the junction becomes, 
in the new geometry, an energy shift equal to the hopping 
t at the atoms in the last stripe of the system: 
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An incoming wave of momentum k and energy efc defines 
a set of coefficients, {a'^}, {Pf}, such that the Ansatz: 
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is a solution of Eq. (P5)) . In the limit A^ ;2> 1, and e^ ^ i, 
these amplitudes are well approximated by Eq. (PT|) . For a 
given energy e/c ~ te~ , we can define only one incoming 



and outgoing set of amplitudes as in Eq. ([27|) . We can 
also define A^ — 1 amplitudes {af }, {/3f } such that 
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(28) 



satisfy Eq. ((25|) . For energies very near the Dirac point, 
the solutions (I27p and ([28]) are such that either the am- 
plitudes a^,af or (3j,f3^ are exponentially small when 

The boundary condition at the edge of the system is 
given by Eq. (j26p where the amplitudes can be expanded 
into an incoming, a reflected, and Mi evanescent waves: 
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The insertion of these expressions into Eq. p6|) leads to A^ 
equations with the A^ — 1 unknowns {C^}, and the reflec- 
tion coefficient R. The even and odd combinations of the 
initial junction problem allows us to define to reflection 
coefficients, R± obtained using the two possible signs in 
Eq. ([26]) . It is easy to show that the transmission coeffi- 
cient of the junction can be written as T = (i?+ — i?_)/2. 

The only solution of Eq.([26|l at low energies requires 
all the amplitudes a? to be exponentially small, as a fi- 
nite set of values a° and (3° are incompatible if Tlj ^ 0. 
Hence, the parameters C^ and R± are determined by 
a set of equations which are independent of the choice 
of sign in Eq. pS)) . so that i?+ w R- with exponential 
accuracy. Hence, the transmission coefficient of the orig- 
inal junction also vanishes with exponential accuracy, in 
agreement with the numerical calculations in Refs.— 1^. 

Results are shown in Fig. [5]. They are consistent with 
those reported in22. The dependence on energy of the 
refiection coefficient shows three characteristic patterns, 
which are repeated as function of the width of the junc- 
tion, as shown in the Figure. This approximate peri- 
odicity is reminiscent of the alternance of metallic and 
non metallic features in carbon nanotubcs and graphene 
nanoribbons. 

The calculation described here requires the existence 
of evanescent waves which have a finite overlap with 
wavepackets derived from both the K and K' points of 
the Brillouin zone, in a similar way to the transmission 
problem analyzed in Ref.^®. The coupling between the 
two valleys can be analyzed in detail by calculating the 
relative weights of the different evanescent waves which 
must be defined at the junction. Results for a junction 
of width A^ = 16 are shown in Fig. [3], where the decay 
length is defined as lev = 1/k in eq. f^B]) . Most of the 
weight is concentrated on evanescent waves with a short 
decay length, which cannot be ascribed to a given valley. 



IV. CONCLUSIONS 

We have studied the transmission, at zero energy, 
through graphene quantum dots attached to leads with 
one incoming and outgoing channels, extending the anal- 
ysis in Refsi^ii. We assume that, for energies at distances 
to the Dirac point smaller than vf/L, where L is the 
typical dimension of the system, extended states can be 
neglected, and the electronic properties are mainly deter- 
mined by evanescent waves, or by localized states induced 
by boundaries. 

We have shown that the dots of various shapes can 
be analyzed, using conformal transformations which pre- 
serve the nature of the states at zero energy. An excep- 
tion is an angle between zigzag boundaries with different 
terminations, where the solution of the problem requires 
the analysis of evanescent waves which do not have a 
well defined valley index, and cannot be described by the 
continuum Dirac equation. 

We find that: 



(i) The transmission vanishes when the leads are at- 
tached to the same edge. 

(ii) The transmission depends exponentially on the ra- 
tio between the size of the dot and the width of the con- 
tact. 

(iii) The shape of the dot changes significantly the 
transmission. 

(iv) The effects of a magnetic field are strongly depen- 
dent on the shape of the dot. A rectangular dot is unique 
in that a magnetic field does not change the transport 
properties (see also^) . In general, the conduction of the 
dot is very sensitive to the magnetic flux through it. A 
similar behavior for dots in the diffusive regime has been 
reported in^S. 

(v) The transmission when the two contacts are at- 
tached to an edge such that both the transmitted and 



reflected wave belong to the opposite valley as the in- 
coming wave depends on the lattice structure at short 
distances. 
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